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ABSTRACT 



Systems which are composed of two or more phases, or subsystems, 
arranged in logical sequence are found frequently in industry and 
defense. Standard procedures for computing lower confidence limits 
on reliability of such systems rely on the use of system data^ 
Engineering changes to any of these subsystems can effect the invali- 
dation of all existing data, necessitating additional, sometimes 
extensive, testing. Such changes are not infrequent in complex sys- 
tems. A need exists for a method of computing lower confidence limits 
on reliability which uses phase data. Some approximation techniques 
have become available. One such technique is currently being used 
by Applied Physics Laboratory, The Johns Hopkins University. Computer 
simulation techniques are used to analyze the accuracy of this pro- 
cedure . 
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CHAPTER I 



INTRODUCTION 

Current technology has produced a class of systems which are com- 
posed of two or more phases, or subsystems, arranged in logical sequence, 
such that the success of the system depends upon the success of the 
component phases. Thus the reliability of the system, R^, depends upon 
the reliability of each of the k phases of which it is composed. We 
express this mathematically as follows: 

k 

R = n R. (1) 

^ i=l ^ 

th 

where R^ is the reliability of the i phase. Standard statistical 
procedures used to compute lower confidence limits (LCL) on R^ 
require the use of system test data. Whenever high degrees of reli- 
ability, at high confidence levels, are demanded of a system, extensive 
testing is required, frequently at considerable cost. 

Engineering changes to any one of the phases may cause the invali- 
dation of all previous system test data. There is a need for a pro- 
cedure to construct lower confidence limits on system reliability using 
phase test data which will salvage as much of the data as possible. 

Approximation procedures are becoming available for so constructing 
these LCL*s, one of which is currently being used by Applied Physics 
Laboratory, The Johns Hopkins University. This procedure is discussed 
in [l].l 

The purpose of this paper is to analyze the accuracy of this pro- 
cedure. Computer simulation techniques will be used to accomplish 
this analysis. 

^Numbers in brackets refer to items in the bibliography with the 
same number. 
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CHAPTER II 



COMPUTATIONAL PROCEDURE EMPLOYED BY APPLIED PHYSICS LABORATORY 

The APL method for constructing LCL*s on R^ assumes the following 

a) mutually independent trials 

b) the distribution of the reliability estimator, R^, is 
normal . 

The formula for computing the 100(l-a)%, one-sided LCL, R^, is 

R^ = R^ - P(l-a) • (2) 

s 

where R^ is the system reliability estimator, P(l-a) is the 

100 (1-a) percentile of the standard normal distribution, and is 

K 

S 

the standard deviation estimator for the R distribution 

s 



k 

R = n R, (3) 

^ i=l " 



where R^ is the phase reliability estimator and is defined by 



R. = s . /n . 
1 11 



(4) 



^ th 

That is, R^ denotes the proportion of successes for the i phase, 

s., out of n. trials. 

1 1 



The value is computed by 

K 



o| = n R. - n 
% i=i ^ i=i 



R. - 

1 



ni - 1 



(5) 



which is an unbiased estimator of the variance of R . Standard devia- 

s 

tion, 5^ , is obtained by extracting the square root of the variance. 
R 

s 

Although biased, this estimator is considered to be sufficiently close 



to OJ 



[ 11 - 



s 

As stated above, one basic assumption of the model is that R^ 
is a normally distributed random variable. As an ancillary result of 
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the simulation, this assumption was verified for those cases for which 
the procedure was reasonably accurate. For each case considered, 1000 
values of were computed and plotted on normal probability paper. 

For one such example see Figure 1. For those cases which proved less 
accurate, the distribution departed from normal. 
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Percentile of System Reliability Estimator 
FIGURE 1 

ONE EXAMPLE OF THE NORMAL DISTRIBUTION OF SYSTEM RELIABILITY ESTIMATOR (CASE 32, TABLE II) 



CHAPTER III 



THEORY OF THE SIMULATION TECHNIQUE 

One standard interpretation which can be applied to the meaning of 

a lower confidence limit involves the method used to compute it. This 

says that regardless of the true value of R^ the computation will 

produce an LCL, R^, which is less than R^ (l-a)% of the time. Use 

can be made of this interpretation in formulating the simulation model. 

If one were to use an exact method of determining R^ and repeated 

the process several times, he could expect that approximately (l-a)% 

of the values would be less than R^. As the number of replications 

approaches infinity, he could expect that exactly (l-a)% of the values 

would be less than R . 

s 

One could use a similar approach to evaluate the accuracy of an 
unproven method. Given a hypothetical system of known R^, and hence 
known R , one could use the proposed method to compute ^ on R . 

If the method was accurate, approximately (l-a)% of the values so com- 
puted would be less than R^ . If the number of replications was suf- 
ficiently large, the percentile of ^ which is less than R^ would 
be a good measure of the accuracy of the procedure. This percentile 

will be denoted A- 

1-a 

The computer simulation technique, used herein, follows the above 

approach. The rapid computational ability of the computer allows one 

to simulate the testing process, compute R^, and replicate the process 

a great many times. One of our measures of effectivess of the method 

will be the relative closeness of A- to R . Other measures of 

1-a s 

effectiveness will be discussed in the results section of this paper. 
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The simulation technique begins by reading from data the parameters 



of the case to be considered; i.e., KK, the number of replications de- 
sired (1000 for all cases considered herein) ; K, the number of phases 
per system; and n^ for all phases; and a. Using a pseudo-random 

number generator and the known phase reliabilities, test data is 
gathered on each of the phases by comparing n^ random numbers with 
(i = 1, 2 , . . • , K) . ^ Each time the random number generator pro- 
duces a number less than or equal to R^ it is counted as a successful 
test* A running total of successes, s^, is kept throughout the course 
of testing. After having so obtained these data, SUBROUTINE RHAT is 
called with parameter a. The values s^ and n^ are made available to 
RHAT via the DIMENSION statement. SUBROUTINE RHAT returns the values 
R^ and R^. This operation is repeated for KK replications to compute 
KK values of R^ and R^ . These values are stored in vector matrices 
for later use. 

SUBROUTINE RHAT is a relatively straightforward application of the 
working formulae contained in Chapter 2. One point worthy of note is 
that only ten values of P(l-a) have been read into computer memory in 
order to conserve memory space. The argument (1-a) is digitized to 
facilitate stowage and retrieval. This limits the values of a which 
can be used to increments of .05, from .05 to and including .50, a 
limitation which can easily be changed as the requirements dictate. 



^The pseudo-random number generator used herein is a library sub- 
routine of the Computer Facility, Naval Postgraduate School, Monterey. 

A printout of this subroutine is included in Appendix II. 

^The reader is invited to follow the program logic by referring to 
Appendix II. Extensive use is made of comment cards to explain the 
operations as they are accomplished. The program notation is intended 
to be the same in meaning as that of the body of the paper. It is, of 
necessity, somewhat different in form because of machine limitations. 

As an aid, a glossary of terms is provided in Appendix I. 
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As the above values of are generated, they are summed for 

eventual division by KK to compute their mean. 

Thus far, the simulation has performed the following operations: 

1) generated the statistics s^ (i = 1, 2 , 3 , . . . , K) 

2) computed the values for the ^ and A1 matrices 

(A1 being an exact duplication of the ^ matrix) 

3) summed the values of R to get the mean, y- 

^ \ 

Next, a sorting routine, SUBROUTINE SHSORT,^ is called for the 

purpose of putting matrices R^ and A1 into ascending numerical 

order. This is done to facilitate plotting these values on probability 

graph paper to determine their distributions. This also facilitates 

obtaining the 100 (1-a) percentile of the ^ distribution, 

Recall the matrix A1 was identical to the matrix before it was 

th 

ordered. Thus, A- equals the KK(l-a) member of the ordered A1 

1-a 

matrix. 

Next, the percentile of which is less than R^ is obtained; 

i.e., the probability ^ is less than R , denoted P(R ). This is 

1j s s 

computed by determining the last element of the ordered A1 matrix 
whose value is less than or equal to R^ and dividing its index number 
by KK. 

The variance of R_ , is calculated next by summing the 

\ L 

squared differences of R_ and y- and dividing by KK, the number of 
replications , 



KK 

- 

1=1 1 



) 



KK 



^This subroutine is a library subroutine of the Computer Facility, 
Naval Postgraduate School, Monterey, Calif. It is printed out in 
Appendix II. 
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The square root is later extracted to obtain the standard deviation, 



The difference between R and A- is then obtained, which is 

s 1-a 

the primary measure of effectiveness of the procedure. 

The remaining statements are editorial in nature; e.g., input/output 
statements which read in data or printout results. Others initialize 
variables, sums, or products to ensure that values are not carried over 
from replication to replication. 

The program output may be summarized as follows: 

a) 1000 values of system reliability estimator, , in 
ascending numerical order 

b) 1000 values of LCL, in chronological order 

c) 1000 values of LCL, R^, in ascending numerical order 

d) system reliability, R^ 

e) 100 (1-a) percentile of the population of R^ population, 

1-a 



f) standard deviation of R^ population, 

g) probability that R^ is less than R^, 

h) mean of ^ 

^ \ 

i) the difference between R and A- 

s 1-a 

j) Zn^q^, a measure of the amount of testing done and the 
unreliability (q^) , or failure rate 
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CHAPTER IV 



SIMULATION RESULTS 

The accuracy of the procedure is measured by the closeness of A^ ^ 

relative to R . A second measure of effectiveness, P(R ), the 

probability that R^ is less than R^ , provides a measure of the true 

confidence level achieved by using this procedure. If the method were 

exact, one could expect A to equal R and P(R ) to equal (1-a) . 

J. oc s s 

In addition to the above criteria, the mean of R^ , , and the 

^ \ 

standard deviation, o- , are also provided as measures of effective- 
ness. These two characteristics provide a measure of the likelihood of 
the procedure producing a lower confidence limit which differs consid- 
erably from R^. It can be shown that for those cases for which the 
procedure is reasonably accurate the distribution of R^ is approxi- 
mately normal. (See Figure 2.) Thus if the procedure is accurate, one 

should expect to find the 1-a percentile to be a known function of 

th 

the mean; e.g., the 90 percentile should be approximately 1.28 stand- 
ard deviations above the sample mean. In addition, the mean and vari- 
ance can be used to construct prediction regions into which R^ is 
likely to fall, assuming again that ^ is normally distributed. 

Several cases were examined in evaluation of this method. In each 
case one or more of the parameters, K, R^, n^, were changed. In 

addition, each case was run for a-values equal to .05, .1, and .2, 
successively o The results of these cases are tabulated in Tables I 
through III. 

The accuracy of reliability testing procedures is most often a 
function of the amount of testing, n^, and the unreliability of the 
subsystems, q^. This quantity is also tabulated in Tables I through 
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Lower Confidence Limit, 
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ONE EXAMPLE OF THE NORMALITY OF THE DISTRIBUTION, R^^ (CASE 32, TABLE II) 



K 

III and is expressed by n.q. . Thus Tables I through III provide 

i=l ^ ^ 

criteria upon which one can judge the accuracy of the procedure and an 
effective measure of when to expect this accuracy. If one knows the 
reliabilities for which he is striving and the amount of testing which 
is economically feasible, he has available a good estimate of the 
accuracy of this procedure. 

The cases considered are tabulated in three separate tables ac- 
cording to the purpose for which they were considered. Cases 1 through 
24, Table I, are the same as those considered by Woods and Borsting [2] 
and are intended to provide a direct comparison with their method. 

Table III contains a ‘'family'’ of cases in which only the parameter n^ 
changes. These cases demonstrate how the accuracy criteria, ^ and 

P(R^), behave as a function of the amount of testing. Figures 3 

criteria 
through 
considered. 



through 8 graphically illustrate the behavior of the accuracy 
as a function of In.q. for the cases considered in Tables I 
III. Each of these graphs utilizes one of the three a-levels 
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MEASURE OF THE ACCURACY CRITERION (R - A^ ) VERSUS Zn. 
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FIGURE 6 

MEASURE OF THE ACCURACY CRITERION [ (1-a) - P(R )] VERSUS En. 
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FIGURE 7 
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MEASURE OF THE ACCURACY CRITERION [(1-a) - P(R )] VERSUS Zn 



SIMULATION RESULTS FOR CASES 1 THROUGH 24 
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in 


VD 


o 


o 


o 


o 


I—H 


CM 


o 


CM 


<3- 


1 


o 


CO 




o 


o 


o 


o 




CJv 


o 


<d- 


00 


r— 1 


o 






o 


o 


o 


o 


<3- 


in 


o 


(3V 


I—H 


< 


• 


CJ^ 


<y\ 


• 


• 


• 


• 




CTv 


• 


00 


CTv 




i-H 


* 


* 


I—H 


1-H 


1-H 


I—H 


* 


• 


irH 


• 


• 




o 






o 






o 






00 








CO 






CM 






CM 












cn 


uO 






I—H 






I—H 






<3- 














m 






m 






O 




















CTi 






CJV 








* 






* 






• 






• 






•H 








o 


















cr 


o 






m 






o 






O 






•H 


Ln 






CM 






«n 






in 








• 






• 






• 






• 






W 


eg 






I—H 






CM 






CM 




























o 








m 






o 






O 






o 








CM 






in 






in 






I—H 


































• 






• 






• 






• 








• 






• 






• 














• 






• 






• 






• 






o 


























o 


CM 




in 


CM 




O 


CM 




in 


CM 






r-l 






CM 


#v 




in 






CM 






•H 




I—H 






1-H 






I—H 






I—H 




a 


II 






II 






II 






II 










II 






II 






II 






II 






•H 






•H 






•H 






•H 








a 


•H 




c 






a 


•H 




a 


•H 


























O 








in 






O 






O 






O 








CM 






m 






m 






I—H 


























0^ 








• 






• 






* 






• 








• 






• 






• 






• 






1 — 1 


• 




I—H 


• 




1 — 1 






I—H 


• 






o 






O 






o 






O 








o 


CM 




O 


CM 




o 


CM 




O 


CNJ 




•H 


« 






• 






• 






• 


0\ 




cr 




I—H 






I—H 






f-H 






1-H 






II 






II 






II 






II 










II 






II 






II 






II 






•H 






•H 












•H 








cr 


•H 




cr 


•H 




cr 


•H 




cr 


•H 






in 


O 


o 


in 


O 


o 


in 


O 


o 


m 


O 


O 




o 


1— 1 


CM 


o 


I—H 


CM 


o 


iH 


CM 


o 


I—H 


CM 


w 


























1 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


S 


& 


O 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


w 


























CO 


























< 


1— H 






CM 






CO 






<r 






o 


eg 






CM 






CM 






CM 
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SIMULATION RESULTS FOR CASES 25 THROUGH 51 





v£) 


vD 


1— 1 


vD 




CM 


CJ^ 




OO 


CM 




00 


o 




m 




00 


vD 


o 


O 




•<3* 


in 


00 


v£) 


v£) 


CM 


r>. 


<3- 


G\ 


!— 1 




I — 1 




m 


r>. 


CM 


in 


<3- 


00 


CJ^ 


O 


O 


<J^ 


oo 




<3- 


vD 


00 


r>. 




a\ 


00 


vD 


00 




in 


00 


vD 


•vf 


r>. 


vD 


D 


O 


o 


O 


O 


O 


o 


O 


o 


O 


o 


O 


O 


o 


O 


O 




r>. 


00 


00 


CJ^ 


•vf 


CM 


m 


oo 




CM 


00 


00 


o 


I — 1 


00 




1— 1 


00 


00 


1 — 1 


t— 1 


00 


00 


r>. 


m 


in 


"d- 




CM 


r>. 


CM 




m 


r>. 


CM 


00 


o 


vD 


r>. 


I — 1 


00 




oo 




O 


00 


00 


vD 


vD 


o 


00 


o 


<r 


in 


m 


a^ 


«— 1 


00 


vD 


<3- 




o^ 




00 




00 




vD 


vD 


00 




r>. 


CJ^ 


00 


00 


a^ 


00 


00 


























O 




o 


O 




o 


O 


o 


O 


O 


O 


o 


o 


O 


o 


o 


o 


o 


o 


o 


w 


o 


CJ^ 


00 


o 


O 


! 1 


CM 


in 


00 




o 


CM 


00 




00 




C3^ 


o 


r— 1 


00 


CM 




CJ^ 


00 


00 


CM 


v£) 




o 


CM 


CM 




'‘d' 


CTk 


00 




a\ 


00 


<r 


00 






00 


r>. 


m 


00 




































S 




a^ 


CM 


00 




vD 


m 




00 


00 




•vf 


i-H 


t-H 


vD 


1 


«— 1 


o 


r>. 


(J^ 


viD 


I— 1 


j— 1 


00 


LO 


00 


CM 


CJ^ 


00 


00 


00 


t— 1 


LO 


t — 1 


CM 


vD 


m 


CTk 


a\ 


v£) 


•vf 




t — 1 


00 


00 


00 


r>. 


< 


o^ 


o 


o 


CM 


o 


O 




O 


o 


00 


CM 


o 


uo 


in 


o 


1 


o 


o 


o 


«— 1 


o 


O 


o 


O 


o 


o 


O 


o 


o 


o 


o 


CO 


































1 






1 






1 


1 


1 


1 


1 


1 


1 


1 


1 




•vf 


1 — 1 


00 


00 


o 


•vf 


CM 


'<!' 


in 




r>. 


t — 1 


o 


o 


m 


a 


in 


m 


vD 


C7^ 




00 


in 


CM 


a\ 


o 


00 


vD 


o 


o 


O 


1 


1 — 1 


m 


00 


m 


00 


C7> 


v£5 


'vf 


t — 1 


o 


00 


I — 1 


o 


o 


C7^ 


t— 1 


v£5 


vD 


vD 


uo 


CM 


t — 1 


in 


vD 


vr> 


o 


00 


CM 


o 


o 


'd* 


< 




00 


00 


00 


r>. 




a^ 


00 


00 


• 




CTn 


• 


• 


C3^ 






















f-H 






I — 1 


1 — 1 






























C3^ 








<t 






o 






oo 






vO 






I — 1 






CO 


vD 






c^ 












CM 






1 — 1 






p^^ 


vD 






CM 






in 






t — 1 






>d- 








00 












00 






a^ 






CJ^ 














' 






* 


















•H 




















o 






o 






cr 


r-. 






o 






in 












00 






•H 


CN 






a^ 






<3- 






vD 












c 


• 






• 






• 






• 






• 








•<r 






00 






•vf 






CM 






r— 1 








t — 1 


vD 


00 


o 


vD 


00 


o 


vD 


00 


O 


vD 


00 


O 


vD 


00 




•vf 


CM 


CM 


>vf 


CM 


CM 




CM 


CM 




CM 


CM 


>d- 


CM 


CM 


•H 
































C 


11 


II 


II 


II 


II 


II 


II 


II 


II 


II 


11 


II 


II 


II 


II 




>— 1 


CM 


00 


I — 1 


CM 


00 


t— 1 


CM 


00 


I— 1 


CM 


00 


t-H 


CM 


CO 




c 


C 


C 


C 


c 


C 


0 


c 


c 


c 


C 


C 


c 


C 


C3 




m 


m 


<3- 








in 


in 


LO 


00 


00 


00 


CM 


CM 


CM 




o 


o 


O 


T 1 


t — 1 


T 1 


o 


o 


o 


o 


o 


o 


O 


O 


o 


































cr 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 




t — 1 


CM 


00 


I 1 


CM 


00 


I — 1 


CM 


00 


1 — 1 


CM 


oo 


t — 1 


CM 


00 




cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 


cr 




m 


o 


o 


m 


o 


o 


in 


o 


o 


m 


o 


o 


in 


o 


o 


o 


1 — 1 


CM 


o 


I— 1 


CM 


o 


I — 1 


CM 


o 


I — 1 


CM 


o 


I-H 


CM 


H 


























• 


• 


• 


i 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


§ 


Q 






s 


a 


25 


25 


25 


25 


25 


25 


a 


a 


a 


a 


H 
































cn 

C 


LTJ 






vD 












00 






CTv 






o 


CN 






CM 






CM 






CM 






CM 
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Table II (continued) 





VO 




CO 


o 




vD 


1—1 






VD 




VD 


o 


CM 


o 


o 




CO 




m 


CO 


CM 


1—1 


CM 


O 




CM 




uo 




1—1 


00 


00 


CO 


uo 


uo 


o 




CNI 


VD 


CM 


o 




O 




1—1 


r— 1 


CJv 


CO 


I— 1 


o> 


rvs. 


VD 


OV 


vD 


C7V 




CO 


VD 


UO 




00 


00 


uo 


VD 


VD 


CO 






CO 


VD 


uo 


CM 


vd- 


CO 


D 


o 


O 


O 


O 


o 


o 


O 


O 


O 


o 


o 


o 


o 


O 


O 


O 


o 


o 




vD 


I— 1 


uo 


CM 


CM 


•<r 


O 


C3V 




o 




CM 


CJ> 


00 


uo 


00 


ov 


o 




O 


CM 


O 




<r 


00 




O 


CM 


CM 




CO 


<r 


00 


00 


VD 


CM 


VD 


J-J 


1—1 


CO 




1 — 1 


CM 


av 


CM 


Ov 


UO 




Ov 


1—1 


Ov 


00 


00 


O 


UO 


O 






CO 






U^ 








1—1 




O 


CO 


Vd- 


ov 


1—1 


uo 


O 


CM 




o> 


ON 


CJv 




VD 


VD 




VD 










ov 


00 


CJV 


CTv 


OV 


OV 


































• 


• 


• 


































o 








O 






O 


O 


o 


o 


O 


O 


O 


O 


O 


o 


o 


o 


o 


o 


O 


CO 




o 


O 


00 




00 


o 


00 




O 






00 


VD 


VD 




o 


m 


PC^ 


1 — 1 


00 


o> 


ov 


CO 


o 


CJv 


•<r 


CJv 


CM 


•<r 


I— 1 


<r 






VD 


CM 


vd- 


' — ^ 




m 


uo 




ov 


00 


•<1- 






m 


Ov 


00 


vd- 




rv. 


vd- 


OV 


rv. 


pL, 






































O 


O 


o 


O 


I— 1 


00 




00 


<1- 


uo 


CM 


1— 1 


CM 


1—1 


I— 1 


I— 1 


1— 1 


<0v 


uo 


1 










CO 


uo 


00 


VD 


■<f 


UO 






o 


o 


O 


O 


00 


CM 


rH 


ov 


ov 


CJV 


o> 


I— 1 


1 — 1 


CM 


00 


CM 


VD 


o 


VD 


CJv 


ov 


C7V 


<Ov 


CO 


00 


<c 


CM 


CM 


CM 


CJv 


o 


o 


CJv 


O 


O 


VD 


I— 1 


O 


vd- 


vd- 


vd- 


vd- 


o 


o 




O 


O 


o 


o 


o 


o 


O 


O 


o 


O 


o 


O 


o 


O 


O 


o 


o 


o 


CO 








































1 


' 


1 


1 






1 




1 


1 






1 


' 


1 


1 




1 










o 


o 






<f 




o 




VD 










o 


vd- 


o 


o 


o 


o 


1— 1 




CM 


VD 


1—1 


CM 


CO 


00 


O 


o 


o 


O 


o 


1—1 


CM 


1 


o 


o 


o 


CO 


CM 


CM 


VD 


UO 


vD 


o 


CM 




o 


o 


o 


o 




OV 


1— 1 


o 


o 


o 








VD 


VD 




<r 


VD 


VD 


o 


o 


o 


o 


vd- 


uo 


<C 


• 


• 


• 


00 






00 






00 






• 


• 


• 


• 


<3v 


OV 




I— 1 


1—1 


1— 1 




















1—1 


I— 1 


I— 1 


I— 1 


• 


• 




o 






00 






00 






00 






CJV 






CJv 








CO 
























<3v 






<Ov 








o 






CO 






CO 






CO 






o 






o 






pci 








! 


















uo 






uo 








Ov 
























Ov 






CTi 














• 






























•H 


O 






O 






O 
























cr* 


O 






O 






uo 






O 






O 






o 






•H 


Ov 






uo 






• 






• 






uo 






uo 






C 


00 






• 






CM 






m 






• 






• 




















1 — 1 






CM 






1—1 






CM 
















uo 






Ln 






uo 






uo 






uo 
































0S 




















• 






• 






• 






0 






• 














• 






• 






• 






• 






• 














• 






• 












• 






• 
























o 




















o 


VD 


CO 


o 


CM 




O 


CM 




o 


CM 




o 


CM 




O 


CM 








CM 


CM 


CO 


•V 




uo 






I— 1 






CO 






uo 






•H 










1—1 






1—1 






1—1 






1 — 1 






1 — 1 




a 


II 


II 


II 


II 






II 






II 






II 






II 
















II 






II 






II 






II 






II 






1—1 


CM 


CO 


•H 






•H 






•H 






•H 






•H 








PJ 


a 


C 


a 


•H 




a 


•H 




0 


•H 




0 


•H 




0 


•H 














UO 






UO 






UO 






m 






uo 


























•V 






0S 




















• 






• 






• 






• 






• 














• 






• 






• 






• 






• 














• 






• 






• 






• 






• 






1 — 1 


1—1 


I— 1 


uo 


4^ 




uo 


0\ 




uo 






I— 1 


#v 




1—1 








o 


o 


o 


o 


CM 




o 


CM 




O 


CM 




O 


CM 




O 


CM 




•H 


• 


• 


• 


• 






• 






• 






• 






• 






cr 










1—1 






1—1 






1—1 






1— 1 






1—1 






II 


II 


II 


II 






II 






II 






II 






II 
















II 






II 






II 






II 






II 








CM 


CO 


•H 






•H 






•H 












•H 








cr 


cr 


cr 


cr 


•H 




cr 


•H 




cr 


•H 




cr 


•H 




cr 


•H 






uo 


o 


o 


in 


O 


o 


uo 


O 


o 


uo 


O 


O 


uo 


O 


o 


uo 


O 


o 


P^i 


o 


I— 1 


CM 


o 


1—1 


CM 


o 


t— 1 


CM 


o 


1 — 1 


CM 


o 


1 — 1 


CM 


o 


I— 1 


CM 


w 






































§ 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


§ 


a 


Q 


S 


3 


<3 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


w 






































cn 






































< 


o 






I— 1 






CM 






CO 






vd- 






uo 






o 


CO 






CO 






CO 






CO 






CO 






CO 
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Table II (continued) 





CM 




CM 


CO 




CM 


VO 




CO 


o 


00 


in 






VD 


cS- 




00 




CN 




1—1 


CO 




CO 


VO 


in 


1 — 1 


00 


o 


m 


in 


<r 


00 




in 


vO 




1— 1 


00 




00 


VO 


00 


O 


o 




'sf 




o 


00 


CM 




CM 


CJV 


CT> 


<p^ 


CM 


CM 


CM 


CM 


in 


<f 


CM 


•vT 


CO 


1—1 


CM 


CM 




00 




VD 


m 


m 


o 


o 


o 


O 


o 


o 


O 


O 


O 


o 


o 


o 


o 


O 


o 


O 


O 


O 


o 




KO 




CM 


r-- 


CO 


00 


o 




<T> 


CTi 


1—1 


o 


CTi 


1—1 


in 


CO 




•‘d- 




in 


o 


<J\ 


tH 


00 


CTi 


00 


00 


VO 


o 


VO 


CM 


o 




VD 


VD 


CO 


in 


J-J 


o 


00 




VO 


CO 


CM 


in 


CO 




vD 


1—1 


O 


VO 


Ov 


'•d' 


VD 


m 


m 


<pd 


m 


1— 1 


CM 






m 




'«d- 


m 




m 


VD 


ov 


VD 


1—1 


CTv 


C7V 


CO 






a^ 


a^ 


(j\ 




cri 


C7V 


<y\ 




CTi 


Ov 


o^ 


in 


"d- 


m 


m 


>d- 


m 






o 


o 


O 




o 


O 








o 














o 




o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


O 


o 


O 


o 


o 


o 


o 


CO 


CM 




VO 




o 


CM 


CO 


CM 


VO 


00 


r>- 


CM 


VD 


CO 


in 


00 


in 


<j\ 


p^i 




\D 


1—1 


1—1 


CO 




in 


'«d- 


CO 


CO 


1—1 


CM 


in 


CO 




CO 


m 




'w' 


in 


00 




m 


m 


m 


CO 












m 


cn 


00 


in 


C7^ 


00 


Pu 






































s 


I— 1 




CTv 


m 


m 


m 


m 


m 


m 


m 


o 








'd- 


VD 


CM 




1 


1— 1 


o 


VD 


















CO 


00 


m 


a^ 


CM 


CM 


1—1 


iH 


CT» 


VO 


1—1 


•*d- 








<r 


•vT 


*<r 


1—1 


"d- 


CO 


(Ti 


o 


CO 




CO 


C 


CM 


o 


O 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


O 


O 


CO 


1—1 


CM 


o 


CM 


1—1 


1 


O 


o 


o 


O 


O 


O 


O 


o 


O 


O 


o 


O 


1—1 


o 


O 


1—1 


o 


o 


CO 








































1 


1 


1 


1 


1 


1 


1 


1 


1 


1 




1 


1 






1 








o 


VD 


00 


o 


o 


o 


o 


o 


o 


o 


in 


m 


o 




CT> 




CM 


in 


« 


1—1 


O 


VO 


o 


o 


o 


o 


o 


o 


o 


m 


VO 


VD 


1—1 




o 


in 


m 


I 


o 


r«- 


CM 


o 


o 


o 


o 


o 


o 


o 


CO 




CM 




r>- 


CM 


'd' 


in 


r— 1 


00 


m 


m 


o 


o 


o 


o 


o 


o 


o 






CO 






o 




00 


c 


a^ 


CT» 


Ov 


• 




• 




• 


• 
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Table II (continued) 
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Table II (continued) 
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SIMULATION RESULTS FOR CASES 52 THROUGH 73 
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Table III (continued) 











iH 




LO 






<3- 


<r 




LO 


i-H 




CM 


00 




CM 




00 


O 




CO 


CM 


CO 


CM 


00 


CM 


CM 


o 


<r 






I-H 


CM 


O 


1^ 


1-1 


rH 


LO 


1— 1 


00 




LO 


as 


OS 


LO 


VO 


os 


1— 1 


o 


rH 


<r 


VO 


VO 




< 


ON 


'sf 


CM 


00 


CO 


1— 1 




o 


a^ 


VO 


00 


00 


vD 


00 


r^ 


LO 




VO 


D 


o 


t— 1 


1— 1 


o 


t— 1 


iH 


o 


1— 1 


o 


O 


o 


o 


O 


o 


o 


O 


O 


O 




LO 


00 


1— 1 


CO 


1— 1 


LO 


VO 


OS 


o 


'■d- 


LO 




00 


o 


LO 


LO 


•d- 


CO 




V£) 


VO 


LO 


VD 


1— 1 


vO 


CO 


CM 


1— 1 


00 


VO 


O 


O 


os 


o 


I-H 


00 


LO 




CM 


00 


o 


VO 


<r 


LO 


VO 


o 




00 


rH 


CM 


a^ 


OS 


CO 


VD 


VO 


VO 




VO 


o 




LO 


CM 


VO 


LO 


CO 




LO 


LO 


CJS 


LO 


LO 


a^ 


LO 


VD 


C3S 




00 






00 






00 






00 






00 






00 


1^ 






o 


O 


O 


O 


O 


o 


o 


o 


o 


o 


o 


O 


o 


O 


O 


o 


O 


O 


W 


vD 


I— 1 


!— 1 


VO 


O 


O 


CM 


CM 


00 


VO 


cjs 


LO 


VO 


00 


o 


I-H 


CO 


CO 




O 


os 


o^ 


00 


O 


CT\ 


1— 1 


CM 


r— 1 


CM 


as 


CM 


VO 


VO 


LO 


<r 


a^ 


C3S 




<1- 


00 


00 


CO 


a\ 


VO 


LO 


00 


00 


LO 


00 






00 


00 


<1- 


00 










































Q 


CO 


m 


LO 


CO 


o 


VO 


CO 


00 


LO 


CO 


CO 


00 


os 


CO 


CO 


o 




d- 


1 


VO 


VO 


o\ 


VO 


r-* 




VO 


CO 


LO 


VO 


00 


00 


CM 


I-H 




CO 


o 


vo 


iH 


CNJ 


CM 


CM 


CM 


CO 


VO 


CM 


o 


<f 


CM 


VO 




a^ 


CM 


rH 


VO 


o 


C3S 


< 


<r 


<r 


o 


•<f 


CO 


O 


o^ 


1— 1 


CM 


O 


CO 


O 


o 


O 


I-H 


00 


CM 


O 




1— 1 


1— 1 


o 


1— 1 


O 


O 


o 


o 


o 


1— 1 


o 


o 


rH 


o 


O 


o 


O 


O 


W 






































Pd 


1 


1 




' 




' 


1 






' 




1 


' 


' 




1 


' 


' 




o 


o 


CO 


O 




CO 




LO 


CM 




o 


LO 




o 


LO 


r>- 


r*^ 


t-H 




o 


o 


<r 


o 


VO 


1— 1 


o 




00 


o 


CM 


CM 


VD 


LO 


VO 


VO 


CO 


O 


1 


o 


o 




o 


CO 




o 




CM 


o 


<!■ 


LO 


VO 


o\ 


LO 


CO 






rH 


o 


o 


LO 


o 


CM 


VO 


LO 


VO 


CO 


VO 


CJS 


VO 


VO 


LO 


<r 






vO 


< 


• 


• 


00 


• 


00 


00 


as 


00 


00 


a^ 


00 


00 


a^ 


00 


00 


os 


00 


00 




1— I 


1— 1 


• 


1— 1 






































































CO 






CO 






CO 






CO 






CO 






CO 






W 


























o- 












Pd 


LO 






LO 






LO 






LO 






LO 






LO 








00 






00 






00 






00 






00 






00 














■ 






■ 












• 






• 






•H 


O 






o 
























o 








O 






LO 






o 






LO 






o 






LO 






•H 


1— 1 






CM 






o 












LO 






CM 






a 


• 






• 






• 






• 






• 






• 






w 


CM 






CM 






CO 






CO 












LO 










CO 






CO 






CO 






CO 






CO 






CO 








n 
























0k 






0k 






<1- 


CM 




LO 


CM 




o 


CnJ 




LO 


CM 




o 


CM 




LO 


CM 






T— 1 






rH 






CM 






CM 






CO 






CO 






•H 




1— 1 






t— 1 






rH 






1— 1 






I-H 






I-H 




a 


II 






II 






II 






II 






II 






II 










II 






II 






II 






II 






II 






II 






•H 






•H 






•H 






•H 






•H 






•H 








a 


•H 




a 


•H 




a 


•H 




a 


•H 




a 


•H 




a 


•H 








CO 






CO 






CO 






CO 






CO 






CO 






LO 






LO 






LO 






LO 


0k 




LO 


r* 




LO 








o 


CM 




o 


CM 




o 


CM 




o 


CM 




o 


CM 




o 


CM 




•H 


• 






• 






• 






• 


0k 




• 


0k 




• 


0k 




cr 




1— 1 






1— 1 






1— 1 






1— 1 






I-H 






I-H 






II 






II 






II 






II 






II 






II 










II 






II 






II 






II 






II 






II 












•H 






•H 






•H 






•H 






•H 










•H 




cr 


•H 




cr* 


•H 




cr 


•H 






•H 




cr 


•H 






LO 


O 


o 


LO 


O 


O 


LO 


O 


o 


LO 


O 


o 


LO 


O 


o 


LO 


O 


o 


Pd 


O 


1— 1 


CM 


o 


1— 1 


CM 


O 


1— 1 


CM 


o 


1— 1 


CM 


O 


I-H 


CM 


o 


I-H 


CM 


W 






































i 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


II 


:3 


S 




Q 




s 


s 




s 


Q 


a 


a 


a 


a 


a 


a 


a 


a 


a 


w 












































































<: 








00 






os 






o 






1— 1 






CM 






o 


LO 






LO 






LO 






VO 






VO 






VO 







33 



Table III (continued) 
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Table III (continued) 
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CHAPTER V 



CONCLUSIONS 

Based upon the results of the simulation, which are tabulated in 
Tables I through 111 and illustrated in Figures 3 through 8, the fol- 
lowing general conclusions can be drawn: 

a) The accuracy of the method is a function of and, in 

general, will increase as this value increases. However, the 
graphs of the accuracy criteria take the shape of a sawtooth 
curve. Therefore, there is no guarantee that increasing 

will improve the accuracy of the method. It is generally true, 
however, that as gets large the difference between 

system reliability, R^, and ^ should converge to zero, 
as should the difference between the desired confidence level, 
(1-a) , and P(R^). 

b) At confidence levels of .9 and .8, and for values of Jln^q^ 

which are sufficiently high, the procedure appears to be rea- 
sonably accurate. Obviously, the numerical value of "suffi- 
ciently high" is subjective and depends upon such unknown param- 
eters as user risk, utility, etc. For most purposes, a value 
of 4.0 should be sufficiently high at both of the 

above confidence levels. The procedure appears to be inaccu- 
rate if R^ is very high or small enough to make 

a failure only moderately probable. 

c) Serious doubt exists as to the acceptability of the method for 
a confidence level of .95 or greater. At the former two con- 
fidence levels there was an apparent convergence toward zero 
error. No such convergence was apparent at the .95 level. 
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Even in Cases 6 and 8, where En.q. was 22,0 and 23.9, 
respectively, the difference, is an order of 

magnitude greater than the corresponding values for a = .1 
and a = ,2. Similarly, the actual confidence level is about 
half of that intended. 



- 37 - 



BIBLIOGRAPHY 



[1] Weiner, S. "Reliability of Multi-Phased Systems," in APL Space 

and Research Programs ^ Status Report #SR lB-2, pp. 83-88. 

[2] Woods, W. Max, and Borsting, J. R. "A Method for Computing Lower 

Confidence Limits on System Reliability Using Component 
Failure Data with Unequal Sample Sizes," Technical Report 
NPS55Wo/Bg8061A, Naval Postgraduate School. 



- 38 - 



APPENDIX I 



GLOSSARY 



Computer Simulation 



Meaning 


Notation 


Variable Name 


System reliability 


R 

s 


RSUBS 


Phase reliability 


R. 

1 


RSUBI(I) 


System reliability estimator 


R 

s 


RHATS 


Phase reliability estimator 


R. 

1 


RHATI 


th 

Successful tests of the i phase 


s , 
1 


s(i) 


t h 

Number of tests of the i phase 


n. 

1 


AN(I) 


Estimator for the variance of the 
reliability estimators 


“r2 

s 


SIGPSQ 


Estimator for the standard deviation 
of reliability estimators 


h 

s 


SIG 


Computed lower confidence limit 


\ 


RHATL 


th 

The (1-a) percentile of the 

computed 


^l-a 


AlA 


Probability that ^ is less than 


P(Rs) 


RSPCN 


th 

The (i) percentile of the standard 

normal distribution 


P(i) 


P(I) 


Mean of the LCL’s generated 


\ 


MEAN 


Variance of the LCL’s generated 




VAR 


Standard deviation of the LCL’s 
generated 


Xj 


SIGMA 


Number of replications per case 


KK 


KK 


Number of phases per case 


K 


K 
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APPENDIX II 



PROGRAM PRINTOUT 



PROGRAM MAIN 



PROGRAM MAIN DOES ALL COMPUTATIONS EXCEPT 
THE COMPUTATION OF THE VALUES RHATS, SYSTEM 
RELIABILITY ESTIMATOR, AND RHATL, LOWER 
CONFIDENCE LIMITS. THESE LATTER COMPUTA- 
TIONS ARE DONE IN SUBROUTINE RHAT 



HERE SET INDEX OF RANDOM NUMBER GENERATOR 

XXX*RNG(C ) 

XXX»RNG(l5n 



HERE SET ALL SUBSYSTEM PARAMETERS AND TEST 
STATISTICS TO ZERO TO AVOID CARRYING VALUES 
OVER FROM ONE CASE STUDY TO THE NEXT. 

DO 9 j»i,ior 
P( J)*G. 

AN(J) =0. 

S(JI»0. 

9 CONTINUE 



HERE READ INTO MEMORY: NUMBER OF REPLI- 

CATIONS DESIRED; NUMBER OF SUBSYSTEMS CON- 
SIDERED; ALFA LEVEL; RELIABILITY AND NUM- 
BER OF TESTS FOR EACH SUBSYSTEM. 

PEA0(5,900)KK,K,ALFA, ( R(J),AN(J), J=1,K) 

9C0 FORMAT (2Il.r, F10.5,/( 2F10.5M 



HERE READ INTO MEMORY TEN VALUES OF STAND- 
ARD DEVIATION OF (1-ALFA) PERCENTILE OF t 
STANDARD NORMAL DISTRIBUTION. 

Pm*l.6A5 

PI2I*1.2016 

P<3»«1.O30 

PI4I-.0A16 

PI5I*.675 

PI6)*.525 

F(7»«.305 

PI0I-.253 

P(9»«.124 

F ( 10»».0 

BK«KK 



HERE INITILIZE THE VALUES OF ALL CUMULA- 
TIVE SUMS AND PRODUCTS TO ZERO AND ONE, 
RESPECTIVELY. 

RSUBS-l. 

00 20 J-l,KK 
RHATLUI «C. 

2C CONTINUE 
PMU» 0. 

SUM- 0. 

SUHA « 0. 
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no nooooooooooooo ooooooooo oooooo ooooo oooooon ooo 



50 



HERE compute SYSTEM RELIABILITY, THAT IS 
THE PRODUCT OF ALL SUBSYSTEM RELIABILITIES. 

rn 5C 17 = 1 , K 
RSUBS=RSUBS»R( IZ) 

CONTINUE 



HERE COMPUTE THE SUM OF PRODUCTS 
( N( n X 0(1) ), THE AMOUNT OF TESTING 
TIMES THF subsystem UNRELIABILITY, FOR 
ALL SUBSYSTEMS. 



DO 45C JR = l ,K 
SUMA = ( AN( JR)>^ (1.- 
450 CONTINUE 



R( JR ) ) )+SUMA 



HERE BEGIN GATHERING 
OF KK REPLICATIONS. 

DO 100 IX=1,KK 



TEST DATA FOR EACH 



in 



HERE SET THE VALUE S(4), THE NUMBER 
SUCCESSFUL TRIALS, TO ZERO FOP EACH 
SYSTEM. BEGIN TESTING. 

CO 10 J=1,K 

5( J)=0. 

CONTINUE 
DO 20C IA=1,K 
NIA=AN( lA) 

DO 2C0 18=1, NIA 



OF 

SUB- 



HERE COMPARE AN(IA) RANDOM NUMBERS TO 
RELIABILITY OF THE (lA)TH SUBSYSTEM, 
AN(IA) BEING THE NUMBER OF TESTS DESIRED 
FOR THE (lA)TH SUBSYSTEM. IF THE RANDOM 
NUMBER IS EQUAL TO OR LESS THAN SUBSYSTEM 
RELIABILITY A SUCCESSFUL TEST IS COUNTED. 



IF ( R(IA)-RNG(1 

199 S( IA)=S( IA)+1. 

200 CONTINUE 



)) 200,199,199 



HERE CALL SUBROUTINE 'RHAT* WHICH COMPUTES 
THE LOWER CONFIDENCE LIMIT FOR THE TEST 
STATISTICS GATHERED ABOVE. THE LOWER CON- 
FIDENCE LIMIT IS RETURNED AND STORED IN 
TWO DIFFERENT VECTOR MATRICES • RHATL • AND 
•Al* FOR FURTHER COMPUTATIONS. THEY ARE 
ALSO SUMMED FOR LATER COMPUTATION OF THEIR 
SUBROUTINE »RHAT« ALSO RETURNS KK VALUES 
OF THE SYSTEM RELIABILITY ESTIMATOR 'RHATS 
THESE values ARE ALSO STORED IN A VECTOR 
MATRIX FOR LATER PRINTOUT. 

CALL RHAT(ALPA,PHIH,CAPA) 

PHATL(IX)=PHIH 
PHATS( IX)=CAOA 
Al(IX)=RHATL( IX) 

PMU=PMU+RHATL (IX) 
icn CONTINUE 



MEAN 



-Al- 



noon ooon ooooo nnnooo noooonnnn nooono ooooo 



39e 

399 

400 



1023 

1024 
1010 



HERE THE LOWER CONFIDENCE LIMIT MATRIX AND 
THE SYSTEM RELIABILITY ESTIMATOR MATRIX 
ARE SORTED INTO ASCENDING ORDER OF THEIR 
VALUES, 

CALL SHS0RT(A1 ,A2,KK) 

CALL SHS0RT(RHATS»CAP1,KK) 

HERE DIVIDE ACCUMULATED SUMS 'OF LOWER CON- 
FIDENCE LIMITS BY THE NUMBER OF REPLICA- 
TIONS TO COMPUTE THEIR MEAN. 

MEAN=PMU/BK 

HERE COMPUTE THE VALUE *KALFA* THAT IS 
( 1 - ALFA ) X THE NUMBER OF REPLICATIONS. 
THIS VALUE IS USED AS THE INDEX NUMBER OF 
THE ORDERED LOWER CONFIDENCE LIMIT MATRIX. 
THUS AKKALFA) IS THE (1-ALFAITH PERCENTILE 
OF THE LOWER CONFIDENCE LIMIT DISTRIBUTION. 

AALFA = ( l.-ALFAM BK 
KALFA=AALFA 
A1A=A1 (KALFA) 



HERE COMPUTE THE 
R SUB S IN A WAY 
DESCRIBED ABOVE. 



PERCENTILE 
SIMILAR TO 



ASSOCIATED WITH 
THE method 



PST=PSUBS 
DO 40C IC=1,KK 

if(rst-ai (icn 

KIC=IC-1 
FIC=KIC 
RSPCN =PIC/BK 
RST=1.1 



398,390,399 



HERE COMPUTE VARIANCE AND STANDARD DEVI- 
ATION OF THE LOWER CONFIDENCE LIMIT DISTRIBUTION. 

NUM=RHATL (ICI-MEAN 

NUMS0=NUM4NUM 

SUM=SUM^NUMSO 

CONTINUE 

VAR = SUM/BK 

SIGMA = SORT(VARI 



HERE COMPUTE THE VALUE RSUBS - A SUB(l-ALFA). 
DIFF = RSUBS-AIA 



HERE 



WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

FORMA 

^F15.5 

FORMA 

FORMA 



6 

[6 

6 

6 

6 

6 

6 

6 

6 

( 

( 

( 



PRINT OUT ALL DESIRED VALUES AND RESULTS. 
,101C ) 

,1011 M J,AN(J» , R(JI,J=1,KI 

,1002IRSU8S 

,10C3IAIA 

,10r4|SIGMA 

,1005IRSPCN 

,1006IMEAN 

,1023IDIFF 

,1C24)SUMA 

/////10X,'R SUB S - A SUB(l-ALFA) EQUALS*, 

/////lOX, 'SIGMA (N(n X 0(11) EQUALS' , FI 5. 5) 
1H1,38X, 'I ' ,12X, 'N( I )' , 9X, 'R( I) '//) 
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1011 FORMAT(38X, I2,F15.0,F15.3) 

1001 FORMAT (/////(10X.5F15.5) I 

1C02 FORMAT (/////10X,'R SUB S (SYSTEM REL lA B I L I TY I • , 

*• EQUALS* fF15. 5) 

1003 FORMAT (///// lOX ,•( 1-ALFA ) X 100 PERCENTILE •, 

**(A SUB 1 - ALFAJ EQUALS ',F15.5) 

ICO^ FORMAT (/////lOX, ♦STANDARD DEVIATION OF COMPUTED ♦, 
♦•LOWER CONFIDENCE LIMIT EQUALS • ,F15. 5 > 

1005 FORMAT ( / // / / 10 X , • P ERCENT I LE ASSOCIATED WITH • 

♦•R sue S EQUALS • tFl5. 51 

1006 FORMAT! ////lOX, 'MEAN OF COMPUTED LOWER CONFIDENCE’, 
• 'LIMIT EQUALS ♦,F15,5) 

1020 FORMAT! ICX. 'lOOC VALUES OF THE LOWER CONFIDENCE ♦, 
♦•LIMIT, LISTED IN ACCENDING ORDER* I 

1021 FORMAThoX, ’1000 VALUES OF THE LOWER CONFIDENCE •, 

♦ 'LIMIT, LISTED IN THE ORDER GENERATED' I 

1022 FORMAT! ItX, '1000 VALUES OF THE ESTIMATOR FOR THE', 
♦* SYSTEM RELIABILITY (RHATS)') 

WRITE!6,102C ) 

WRITE!6,10C1 )A1 
WRITE!6,1C2U 
WRITE!6,10C1)RHATL 
WRITE!6,1022) 

WRITE!6,:001 )RHATS 

STOP 

END 
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SUBROUTINE RHAT 



SUBROUTINE RHAT COMPUTES THE VALUES RHATS 
(SYSTEM RELIABILITY ESTIMATOR) AND RHATL 
(LOWER CONFIDENCE LIMIT) USING AS ITS 
CALLING PARAMETERS THE STATISTICS GENERATED 
IN PROGRAM MAIN. IT IS A STRAIGHTFORWARD 
APPLICATION OF THE APPLIED PHYSICS LAB- 



COMMON R(IOC) f AN( 100) ,S(100)fPHATL(1000),F(l0), 
♦Al (1000) ,A2 (1000) ,KK,K 



HERE INITIALIZE VALUES AS NECESSARY 

PIPS0=1. 

PR0D»1. 

CAPP=1, 



HERE 'RHATP IS SYBSYSTEM RELIABILITY ESTIMATOR 

DO 50C JA=1,K 
RHATI=S(JA)/AN( JA) 

FSO*RHATI»RHATI 



HERE PIPSO IS THE PRODUCT OF THE SQUARED 
SUBSYSTEM ESTIMATORS. 

PIPSO=PIPSO*PSO 



THE PURPOSE OF THIS STATEMENT IS TO PRE- 
VENT EXPONENT UNDERRUN, THAT IS, IT PRE- 
VENTS NUMBERS FROM GETTING SMALLER THAN 
THE COMPUTER CAN EFFECTIVELY HANDLE. IT 
SETS VERY small NUMBERS EQUAL TO ZERO. 

IF( PIPSQ-.OOCOOCC 1 ) SOI ,5C2,502 
501 PIPSO=0. 

PROD=C. 



THIS STATEMENT IS A CODIFICATION OF THE 
SECOND TERM OF THE RIGHT HAND SIDF OF 
EQUATION (5), CHAPTER II. 

502 PROO=PROD»(PSQ-(RHATI* (l.-RHATI) ) / ( AN ( J A )- 1. ) ) 

CAPP*CAPP*RHATI 
500 CONTINUE 



THIS STATEMENT IS A CODIFICATION OF THE 
ENTIRE EQUATION (5), CHAPTER II. 

SIGPSO=PIPSO-PROD 

SIGP*SORT(SIGPSQ) 

BALFA»ALF *20. 

JALFA»8ALFA 



HERE COMPUTE THE VALUE OF THE LOWER CONFI- 
DENCE LIMIT, CALLED PHA IN THIS SUBROUTINE. 

PHA *CAPP-(P( JALFA)*S IGP) 

RETURN 
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SUBROUTINE SHSORT 



SUBROUTINE SHSORT ( A ,KEY ,/ N/ » 
DIMENSION A(NI,KEY(N» 

Ml«l 

6 M1«M1*2 

IF<M1-N» 6,6,8 
8 Ml*Ml/2-l 

MM«MAX0(Ml/2,l) 

GO TO 21 
2C MM»MM/2 

IF(MM)100,100 ,21 

21 K»N-MM 

22 DO 1 J*1,K 
II«J 

11 IM=II+MM 

if(a(im»-a( 1 1 n 30,1,1 
30 TEMP=A( 1 1 » 

IT»KEY( in 
A< I n=A{ IM» 

KEY( in=KEY( IMI 
A( IMI=TEMP 
KEY( IM)=IT 
II*I I-MM 

iF(in 1,1,11 
1 CONTINUE 
GO TO 20 
100 RETURN 
END 



- 45 - 



FUNCTION SUBROIJTINE RNG 



FUNCTION RNGIN) 

NR=N 

IF(NRI10,lCt20 
10 IX=30517 
NR=NR+l 

20 DO 50 1=1, NR 
IY=IX*65539 
IF(IY)5.6,6 

5 IY=IY+2147A836<f7^1 

6 RNG=IY 

RNG=RNGY.4656613E-9 
50 IX=IY 
RETURN 
END 
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